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ABSTRACT 


The  one-dimemionai  inui  flow  through  s typical  one-cell  built-up  bcx 
beam  was  investigated  experimentally  to  determine  the  stresses  and  distortions  which 
result  from  a transient,  non-uniform  temperature  distribution.  The  model  was  exposed 
to  u radiant  heater  which  produced  the  desired  maximum  temperature  of  400°F  in  the 
nearest  surface  within  two  minutes.  Eleven  temperature  and  nine  strain  gags  readings 
were  recorded  continuously  during  the  interval.  A comparison  of  the  Mins  historic:  of 
the  temperature  distribution  with  the  bax  beam  interior  empty  and  then  with  it  filled 
with  insulation  show  that  conduction  is  the  primary  mode  of  heat  transfer.  Graphical 
and  tabular  comparisons  between  rh±  experimentally  determined  stresses  and  the  stresses 
calculated  on  the  basis  of  the  temperatures  measured  at  eleven  points  are  presented  for 
heating  times  of  30,  60,  90  and  120  seconds.  The  agreement  is  fair  for  early  heating 
times  and  poor  for  later  times.  Conditions  at  the  lower  flanges  and  thermal  buckling 
may  explain  the  discrepancies. 
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INTRODUCTION 


The  development  of  the  atomic  bomb  and  the  advent  of  supersonic  flight 
have  raised  the  phenomenon  of  heat  conduction  to  a very  high  position  In  the  list  of 
problems  faced  by  aeronautical  designers.  Security  measures  have  restricted  the  amount 
of  Informption  on  the  effects  of  thermal  radiation  from  atomic  bombs  on  aircraft  that  Is 
generally  available.  But  a sufficient  amount  can  be  gleaned  from  various  sources  (see 
Refs.  1 and  2)  to  realize  that  the  calculation  of  temperature  distributions  and  thermal 
stresses  may  be  necessary.  An  excellent  series  of  papers  (Refs.  3,  4 and  5)  have  ana- 
lyzed the  time  dependence  of  the  temperature  distribution  In  thin  solid  wings  at  high 
Mach  numbers.  The  NACA  has  reported  analytical  and  experimental  investigations 
cf  the  e?£»r-»*  of  non-uniform  temperature,  distributions  on  semi-monocoque  structures 
and  their  components  (Refs.  6,  7,  8 and  9).  Reference  iC  present:  a st'.Hy  of 
aspect  of  thermal  stresses  In  aircraft  structures  due  to  aerodynamic  heating.  The  above 
ten  references  represent  only  a portion  of  the  pertinent  ones  which  are  available  in  this 
field.  Comparatively  few  published  reports,  however,  contain  experimental  data  on 
the  temperature  and  stress  distributions  in  any  type  of  aircraft  structure  which  Is  sub- 
jected to  a rapid  rate  of  heat  input.  It  Is  the  principal  purpose  of  the  report  to  des- 
cribe and  report  the  results  from  some  experiments  of  this  nature. 
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DESCRIPTION  OP  EXPERIMENTAL  METHOD 


The  purpose  of  this  Investigation  is  to  study  experimental  ly  the  manner  in 
which  heof  is  conducted  through  a one-cell  box  beam  and  to  determine  the  stresses  and 
distortions  which  accompany  the  resulting  transient,  non-uniform  temperature.  By  de- 
signing the  box  beam  with  a fairly  large  span-to-depth  ratio,  the  local  effects  of  the 
free  end  and  of  the  root  were  made  negligible.  As  a result,  the  heat  flow  was  essentially 
a one-dimensional  phenomenon,  and  the  measured  stresse  s could  then  be  compared  with 
the  values  obtained  from  a simple  analysis. 


2. 1 The  Model 


The  experimental  model  in  this  test  program  consisted  of  a one-cell  box 
bourn,  which  was  3 inches  deep,  6 inches  wide  and  25  3/4  inches  long.  The  box  was 
belted  together  with  24  S-T  Phillips  flat  heed  machine  screws  and  24  S~T  hexagonal 
nuts.  An  inch  and  three  quarters  of  the  length  at  one  end  of  the  model  was  used  to 
mount  heavy  steel  angle  irons  that  fastened  the  model  to  the  bed  plates  of  an  erector 
ser  In  a cantilevered  position  (see  Figure  2.  1).  Although  the  cross  section  of  this  box 
beam  is  simple,  several  problems  hod  to  be  solved  before  the  design  assumed  the  form 
shown  in  Figures  2. 2,  2. 3,  2. 4 and  2. 5.  The  design  had  two  basic  aims:  (1)  to  dup- 
licate a typical  one-cell  built-up  box  beam  consisting  of  spar  caps,  shear  webs  and 
skin  and  (2)  to  minimize  the  number  of  contact  surfaces  without  undue  fabrication  effort 
in  order  to  gain  better  heat  transfer  between  the  web  and  the  spar  cap.  The  shear  webs 
and  the  soar  caos  were  simulated  by  an  0.G65  inch  thick  24  S-T  aluminum  alloy  sheet 
and  a 3 inch  61  S-T  structural  alui.Inum  channel,  respectively.  Aluminum  sheet  of 
0.065  Inch  thickness  wets  used  because  of  its  availability  and  also  because  the  desired 
duration  of  the  heat  input  could  not  be  obtained  without  seriously  damaging  the  thinner 
skin.  A structural  aluminum  channel  was  used  because  it  was  easily  obtained  from  cur- 
rent stock  and  because  the  mechanical  properties  of  61  S-T  are  not  very  different  from 
those  of  24  S-T  aluminum  sheet  (see  Table  A.  1 in  Appendix  A).  The  cross  section  of 
the  channel  is  shown  In  Figure  2.5.  The  shaded  portions  indicate  the  material  which 
was  removed  in  order  to  simulate  a perfect  ioint  between  the  shear  web  and  spar  cap. 

The  one-inch  spacing  of  the  bolts  which  corrected  the  skin  panels  to  the 
channel  was  chosen  sc  as  to  preclude  the  possibility  of  inter-bolt  thermal  buckling. 

The  stresses  in  the  skin  panels  were  calculated  on  the  basis  of  an  assumed  cubic  tem- 
perature variation  across  the  depth  of  the  beam  with  an  assumed  maximum  temperature 
for  the  bottom  skin  of  400°  F.  These  calculated  stresses  were  then  used  to  compute  the 
spacing  necessary  to  prevent  buckling  between  the  bolts. 
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PIGURE  2. 1 EXPERIMENTAL  SETUP 


FIGURE  2.2  EXPLODED  VIEW  OF  MODEL 
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During  the  test  program,  tha  model  was  held  in  a cantilevered  position  by 
means  of  a pair  of  heavy  steel  angle  irons  which  were  connected  to  the  outer  skin  on 
each  side  of  the  model  by  3/8  inch  steel  bolts.  The  bolts  extended  through  the  angle 
Irons  and  the  skin  of  the  model  and  into  a one  Inch  thick  steel  block  which  was  fitted 
to  the  inner  dimensions  of  the  model.  The  other  legs  of  the  angle  irons  were  fastened 
to  heavy  bed  plates  which  were  formed  from  an  erector  set  in  the  laboratory.  This  rigid 
support  for  the  box  beam  eliminated  cny  chance  of  unwanted  displacements  of  the  mods! 
(see  Figures  2. 1 and  2. 3). 


2.2  Instrumentation 

The  temperature  distribution  in  the  box  beam  was  fcu:»  by  using  eleven 
Ruge-DeForest  RDF  Stlkons,  Type  BN-3  resistance  thermometer  elements,  mounted  at 
the  midspan  station  as  shown  In  Figure  2.6  and  designated  by  the  symbols,  T-i  to  T-i  i. 
These  gages  were  found  to  be  well  suited  to  transient  measurements  because  of  their  low 
thermal  inertia.  Thetr  sensitivity  to  strain  is  negligible,  and  their  operating  range  ex- 
tends up  to  400°F,  a I though  at  higher  temperatures  the  gages  could  be  used  oniy  for 
short  periods  of  time. 

The  strain  pickups  used  on  the  tower  skin  and  shear  web  were  Baldwin  SR  -4 
temperature  compensated  (from  + 50°  to  + 250°F)  strain  gages,  type  EBDF-7D  designed 
specifically  far  the  use  on  aluminum  alloy.  It  should  be  noted  that  the  type  EBDF-7D 
gages  give  outputs  which  are  proportional  to  the  quantity, 

€-c?zrr 


where 


€ is  the  total  strain 

(X  is  the  coefficient  of  thermal  expansion 

»*  the  change  In  temperature  from  its  equilibrium  value. 

The  above  quantity  represents  only  the  strain  due  to  stress  ard,  lienee,  the  gages  do  not 
read  strain  in  the  conventional  sense  of  the  word.  The  ccefrciont  of  thermal  expansion 
must,  ef  course,  be  the  samu  os  that  for  the  material  under  test.  The  upper  skin  had 
only  Baldwin  SR-4  bake  life  base  strain  gages.  Type  AB-7,  because  ot  the  small  tem- 
perature changes  expected.  It  was  founa,  however,  that  even  though  the  temperature 
changes  in  the  AB-7  gages  were  small,  a temperature  correction  had  to  be  applied  to 
their  strain  readings  because  the  dummy  gogei  used  In  the  bridge  circuit  were  wire- 
wound  resistors  which  operated  at  room  Temperature.  No  correction  fur  gag®  factor  due 
to  ramperuhn*  changes  was  nscetsory  for  either  type  of  strain  gage.  The  EBDF-7D  gages 
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ere  tenrperature  conpensated,  while  the  gage  factor  of  the  AB-7  is  constant  ft*  tem- 
peratures up  to  300°F  (Ref,  20).  All  of  the  w.hove  ;/3ges  were  installed  according  to 
the  manufacturer's  recommended  procedures  and  baking  cycles.  The  strain  gages  are 
designated  by  the  symbols,  S-l  to  S-l  1,  and  their  locations  are  shown  in  Figure  2.6. 

A twelve -channel  Consolidated  Engineering  Company  recording  oscillo- 
graph was  used  to  record  the  temperature  changes  at  the  slowest  possible  film  speed  of 
1/8  of  an  Inch  per  second  so  as  to  minimize  the  amount  of  film  that  had  to  be  developed. 

A special  set  of  dummy  resistors,  used  in  conjunction  with  the  bridge  circuits  of  the  oscil- 
lograph, reduced  the  amplitude  of  the  signal  received  from  the  Stikon  gages.  Bald- 
win Portable  Strain  indicators.  Type  L,  measured  the  strain  gage  reading  directly  in  micro* 
inches  per  inch.  Figures  2.7  and  2.  8 present  block  diagram  sketches  of  the  elements  of 
each  circuit.  Figure  2. 1 0 shows  the  test  equipment  which  was  used.  In  addition,  two 
dial  gages,  shown  in  Figure  2, 1,  were  mounted  at  the  free  end  of  the  model  for  the  pur- 
pose of  measuring  deflections.  These  gages  were  used  to  mark  the  beginning  of  the 
nouiif^  cycle  and  — :;c  Indicated  total  d*fUet5on  at  the  end  of  the  heating  period. 


2. 3 The  Radiant  Furnace 

An  investigation  of  both  commercial  and  special  purpose  heaters  disclosed 
that  the  heater  of  Reference  2 would  be  adequate.  This  radiant  furnace  originally  con- 
sisted of  twelve  heating  alements  of  open  coiled  windings  cf  nichrome  wire  which  were 
mounted  in  an  insulated  aluminum  reflecting  box,  6 Inches  wide  and  24  inches  long. 

Each  element  is  composed  of  0. 064  inch  diameter  wire  and  has  a resistance  of  36  ohms. 

The  heating  coils  were  connected  in  parallel  and  nvsunted  on  heavy  buss  bars  located 
at  the  ends  of  the  reflecting  box.  Additional  supports,  consisting  of  transite  separators 
spaced  approximately  3 1/2  inches  apart,  proved  necessary  to  minimize  the  sag  of  the 
coi  Is. 

An  87  KW  D.C.  generator,  which  normally  drives  the  M. I.T.  Flutter  Tunnel 
also  usually  supplies  the  power  required  by  the  heater  (see  Figure  2. 9).  But  since  the 
tunnel  was  in  operation  during  most  of  this  investigation,  the  heater  had  to  be  modified 
so  as  to  operate  on  the  electric  power  provided  by  *he  Cambridge  Electric  Power  Com- 
pany. The  alternate  power  available  in  the  Structu  is  Laboratory  was  supplied  by  a 
220  volt,  100  anpere,  three  phase,  Y -connected  transermer.  Modification  of  the 
heater  consisted  of  separating  the  twelve  heating  coils  into  three  equal  groups  to  ferm 
a Y-type  three-phase  circuit.  This  alteration  v,*as  accomplished  very  easily  by  dividing 
the  buss  bar  at  one  end  into  three  sections  with  four  coils  per  section  (see  Figure  2. 9). 

In  the  original  construction  of  the  heater,  the  resistance  of  each  coii  was  considered 
equal  with’n  0.01  ohm.  Thus,  It  was  assumed  that  there  would  be  no  unbalance  of 
power  in  the  circuit.  Since  the  current  varied  only  0. 1 of  an  ampere  between  branches, 
the  assumption  was  justified.  In  addition,  the  current  in  each  branch  remained  essentially 
constant  (13.5  amperes)  during  the  heating  period. 
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FIGURE  2.7  BLOCK  DIAGRAM  OF  TYPICAL  STRAIN  GAGE  MEASURING  CIRCUIT 
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FIGURE  2,  8 BLOCK  DIAGRAM  OF  TEMPERATURE  RECORDING  CIRCUIT 
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FIGURE  2. 10  RECORDING  EQUIPMENT 


I Previous  tests  showed  the  heater  capable  of  supplying  up  to  20,  COO  watt: 

j (Rff f.  2).  With  above  configuration,  however,  the  hearer  supplied  only  5, 200  watts, 

j The  heater  rested  on  a small  table  and  was  directly  under  the  box  beam  with  an  average 

| distance  of  7/E  of  an  inch  between  the  box  beam  and  the  heating  coils  (see  Figure  2. 1). 

The  power  output  could  have  been  increased  either  by  connecting  more  r.ichrcme  heating 
colls  in  parallel  in  each  branch  of  the  three  phase  circuit,  or  by  connecting  the  heater 
as  a delta,  thrso  phase  circuit.  This  increase  proved  unnecessary,  however,  since  the 
desired  maximum  temperature  of  400°F  in  the  bottom  skin  of  the  modal  was  obtained 
within  two  minutes  of  heating  time  for  the  heater  in  its  present  configuration.  The 
i power  of  5.2  KW,  supplied  to  die  coils,  is  equivalent  to  approximately  5. 0 BTU/sec 

of  radiant  heat  energy  on  the  assumption  that  there  are  no  losses.  The  absorption  by 
the  mode!  depends  upon  a number  of  factors  such  as  absorptivity  of  the  aluminum,  the 
distance  from  coils  to  skin  and  convective  losses.  There  is  reason  to  believe  that  less 
than  50%  of  the  electrical  power  input  was  absorbed  by  the  model,  but  a mare  detailed 
analysis  and  expert  mental  study  of  such  losses  was  deemed  beyond  the  scope  of  this  in- 
vestigation. 


2. 4 The  Test  Procedure 

The  peak  temperature  for  any  portion  of  the  box  beam  was  limited  to  400°F 
because  at  higher  temperatures  the  material  properties  begin  to  decrease  very  rapidly 
(see  Figure  2. 11).  This  limitation  required  the  selection  of  a heating  time  which  would 
not  overheat  any  point  of  the  model  but  which  would  still  be  of  sufficient  duration  to 
induce  appreciable  thsrmaelastic  effects.  Calculations  indicated  that  the  bottom  skin 
would  reach  350°F  in  approximately  one  minute  if  all  the  electrical  energy  o?  the  radiant 
furnace  at  its  rcfecl  output  of  5. 2 kilowatts  were  absorbed  by  the  bottom  skin  alone.  A 
more  realistic  duration  of  two  minutes  was  determined  by  the  following  simple  experiment. 
A fiat  shee*  of  24  S-T  aluminum  alloy  was  placed  over  the  furnace  and  subjected  to  vari- 
ous iRsutlng  dictions.  From  measurements  with  copper  constantin  thermocouples  mounted 
on  the  sheet,  it  was  found  that  a period  of  approximately  two  minutes  was  required  to 
raise  the  temperature  of  this  test  sheet  to  350°F. 

To  define  further  the  capabilities  of  the  furnace,  the  aluminum  test  sheet 
was  heated  continuously  until  a steady-state  condition  was  attained.  This  steady-state 
condition  was  reached  after  8 minutes  of  heating  time  with  the  sheet  at  a final  tempera- 
ture of  650°F.  In  this  experiment  free  convection  took  place,  arid  the  opposite  side  of 
the  sheet  radiated  to  the  rest  of  the  laboratory.  Similar  conditions  do  not  prevail  in  a 
box  beam  since  the  interior  surfaces  are  radiating  to  a relatively  small  volume  and  free 
convection  in  the  interior  is  also  restricted.  The  importance  of  heat  transfer  by  radia- 
tion and  convection  in  the  interior  of  the  box  was  investigated  by  comparing  measured 
temperatures  and  strains  for  two  different  configurations,  namely,  on  empty  interior 
end  an  inferior  filled  with  batts  of  insulation. 


An  actual  test  run  proceeded  Jr;  the  following  monr**?*"-  <Kr 
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to  the  test  equipment  were  turned  on  to  allow  the  equipment  to  warm  up.  The  eleven 
channels  of  the  temperature  recording  equipment  'were  balanced  at  Identical  signal 
levels,  and  the  room  temperature  displacements  were  recorded  for  each  channel.  The 
portable  strain  indicatcrs  were  balanced,  and  their  initial  reading;  were  recorded  since 
these  instruments  are  mill  devices.  Only  four  of  the  nine  strain  gages  on  the  model 
could  he  read  during  each  test  run  because  of  limitations  In  the  equipment  available. 

It  was  found  that  the  temperature  and  strain  distributions  were  similar  for  each  test  run, 
and  the  strain  readings  used  in  the  calculations  are  average  strains  for  at  least  five  tost 
rum  for  each  configuration.  An  asbestos  insulating  blanket  was  placer!  over  the  model 
and  heater  to  minimize  the  effects  of  convection  currents.  Before  each  tesrt  run,  the 
room  temperature  was  read  from  a bulb  thermometer  and  recorded.  It  was  necessary  to 
know  mo  room  temperature  because  it  was  the  Initial  teinperature  of  model  which.  In 
turn,  was  also  the  reference  temperature  for  the  temperature  gages.  As  it  happened, 
the  room  temperature  remained  almost  constant  during  the  test  program.  The  model 
*ve«  insulated  from  the  heater  by  a transite  sheet  placed  between  them  until  the  heater 
reached  the  desired  temperature.  Then  the  transire  was  removec!. 

When  eH  w"?  in  readiness,  the  power  for  the  electric  heater  was  turned  on, 

and  the  dial  gages  attached  to  the  model  were  wutched  for  the  initic!  deflection.  As 
soon  as  they  began  to  deflect,  the  transite  insulating  sheet  between  the  model  and  the 
heater  was  withdrawn,  and  the  timing  period  started.  A stop  watch  was  used  to  indicate 
ten  second  intervals  during  the  two  minute  heating  period;  at  these  times  the  strain  gages 
were  read  directly.  An  internal  timing  pulse  was  used  to  Insert  10  second  reference  marks 
directly  on  the  oscillograph  flim  using  the  twelfth  channel,  while  the  other  eleven  channels 
were  ccnthxiously  recording  the  readings  of  the  eleven  temperature  gages.  At  the  end  of 
the  two  minute  heating  cycle,  the  tip  deflections  of  both  dial  gages  were  recorded.  The 
power  to  the  heater  was  shut  off  and  an  electric  fan  near  the  model  was  turned  on  to  aid 
in  cooling  the  model  to  room  temperature.  This  process  required  a half  an  hour. 

After  the  completion  of  tests  on  the  model  with  and  without  the  insulation, 
the  temperature  recording  equipment  was  calibrated.  Each  of  the  eleven  temperature 
gages  on  the  model  were  replaced  in  the  oscillograph  channel  circuits  by  a decade  box. 
Then,  known  changes  in  resistance  were  applied  through  the  amplifier  circuits,  and  the 
corresponding  deflections  were  recorded  by  the  oscillograph  for  the  same  attenuations  of 
the  amplifiers  that  were  used  during  the  test  runs.  This  procedure  produced  a relation 
between  ths  change  in  galvanometer  deflection  and  change  in  resistance  for  a particular 
attenuation.  The  deflections  of  all  of  the  temperature  gages  during  the  rest  runs  were 
converted  to  resistance  changes,  and  since  the  Initial  temperature  of  the  inode!  was 
known,  the  total  resistance  could  be  determined.  These  total  resistances  were  converted 
to  temperatures  by  using  the  tempered*  gags  calibration  charts  supplied  by  the  manu- 
facturer. 
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THt  EXPERIMENTAL  RESULTS 


The  experimental  results  consist  of  ('he  reading:  obtained  from  eleven  tem~ 
perature  gages  and  nine  strain  gager*  which  were  located  at  tho  mid-span  station  as  shown 
in  Figure  2.6.  The  strain  readings  ’were  recorded  directly  in  microinches  whereas  the 
temperatures,  recorded  i..  analogue  form  on  oscillograph  paper,  Jxxl  to  be  converted 
into  temperatures  by  calibration  charts. 


I 


3. 1 The  Temperature  Distribution 

Figure:  3. 1 and  3, 2 present  the  transient  temperature  curves  for  the  eleven 
temperature  gages  with  uninsulated  and  insulated  box-beam  interiors,  respectively.  The 
same  data  Is  tabulated  in  table  A.  2.  A eonpuniwi  of  the  two  figures  shows  that  the 
temperatures  for  the  uninsulated  Configuration  1 are  approximately  seven  to  twelve  de- 
grees higher  than  those  at  corresponding  points  on  the  insulated  Configuration  2 at  the 
end  of  heating  period.  The  smallest  differences  both  in  absolute  and  percentage  terms 
occurred  along  the  bottom  skin.  On  the  other  hard,  the  top  skin  with  the  smallest  tem- 
perature increases  had  the  largest  differences  both  in  absolute  terms  and  in  terms  of  the 
percentage  of  the  temperature  rise.  Thus,  the  rock  wool  insulation  appears  to  have  two 
effects:  (I)  the  prevention  of  heat  transfer  across  the  interior  whether  by  radiation  or 
convection  and  (2)  the  retention  of  a portion  of  the.  thermal  energy  which  it  receives 
from  the  interior  surfaces  of  the  box  through  conduction  and  radiation,  i.e.,  a heat 
sink.  The  first  effect  would  account  for  the  fact  that  the  temperature  differences  pro- 
duced by  the  insulation  are  larger  in  the  top  skin  than  in  the  bottom  skin,  in  the  un- 
insulated case,  heat  rccche:  the  top  skin  bv  conduction,  convection  and  radiation,  and 
the  addition  of  the  insulation  suppresses  the  latter  two  modes  of  heat  transfer. 

Figures  3.  3 and  3. 4 show  the  time  variation  of  the  ten^peraturs  distribution 
in  the  web  of  the  box  beam  for  Configurations  i and  2,  respectively.  Figure;  4. 16  and 
4. 17  present  composite  curves  of  the  temperature  distribution  In  the  model  for  the  in- 
sulated case  at  four  different  times.  These  curves  are  plotted  in  a peripheral  coordinate 
system  which  originates  at  the  mid-point  of  the  bottom  skin.  The  increments  of  tem- 
perature increase,  , are  plotted  instead  of  the  absolute  values  of  the  temperature, 
and  an  assumed  mode  of  temperature  is  also  plotted  for  purposes  of  comparison  (see  Sec- 
tion IV). 
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FIGURE  3.2  VARIATION  OF  TEMPERATURE  VS  TIME  FOR  CONFIGURATION  2, 
(ROCK  WOOL  INSULATION) 


- 19  - 


Tr-'r— - 


ASRLTR  25-13 


3.2  The  Distribution  of  the  Strain  d us  to  Stress 


Figures  3.5  and  3.6  give  the  experimental  time  variation  of  the  stress  at 
the  nine  strain  gage  locations  for  the  twojeonfigurations.  The  same  data  is  tabulated 
in  Table  A.  3.  The  largest  strains  are  compressive  and  occur  in  the  bottom  or  heated 
skin.  The  largest  tensile  strains  occur  in  the  web. 

The  time  variation  of  the  strain  distribution  in  the  web  is  shown  in  Figure 
3.7.  For  the  first  50  seconds,  the  strain  distribution  is  essentially  linear,  but  after 
that  time  interval,  the  distribution  in  the  web  becomes  increasingly  complex.  If, 
however,  only  a narrow  portion  at  the  center  of  the  web  is  considered,  i.e.,  the 
middle  1. 20  inches,  the  strain  distribution  is  observed  to  progress  from  a linear  dis- 
tribution ta  a gradually  curved  distribution.  The  character  of  the  strains  changes 
rather  drastically  in  the  distance  between  this  central  portion  and  the  flanges.  Zero 
strain,  however,  occurs  at  about  the  same  point  for  ail  times.  Figure  3. 8 presents 
the  time  variation  of  the  strain  in  the  hat'tewn  skin.  If  a linear  distribu- 

tion is  assumed,  the  strains  at  various  times  intersect  at  approximately  the  same  loca- 
tion. 


3.  3 The  Deflection  at  the  Free  End 


The  fact  that  the  readings  of  the  two  dial  gages  never  differed  more  than 
0. 002  inches  affords  evidence  that  only  negligible  twisting  of  the  model  occurred 
during  the  test  runs.  This  agreement  between  the  gages  gives  a fairly  good  indication 
that  the  temperature  distribution  in  the  beam  was  symmetrical  obout  the  vertical  center 
lire,-  although  the  evidence  is  not  conclusive.  Also,  the  maximum  displacement  at 
the  end  of  each  run  was  within  0. 005  inches  of  the  average  for  all  runs.  The  average 
displacement  for  Configuration  1 was  0.  168  inches  and  for  Configuration  2,  0. 161 
inches.  This  bvidsr.es,  together  with  the  close  agreement  of  the  measured  tempera- 
tures, offers  a fairly  good  indication  that  the  test  conditions  were  duplicated  in  the 
series  of  runs  required  by  the  limited  number  of  strain  reading  device*  (see  Section 
2.4). 
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FIGURE  3.5  EXPERIMENTAL  TIME  HISTORY  OF  STRESS  AT  STRAIN  GAGE  LOCA 
TIONS,  CONFIGURATION  1 
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FIGURE  3.6  EXPERIMENTAL  TIME  HISTORY  OF  STRESS  AT  STRAIN  GAGE  LOCA- 
TIONS, CONFIGURATION  2 
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ANALYSIS  OF  RESULTS 


in  principle,  the  calculation  of  the  three-dimensional  distribution  of  ther- 
mal stress  Is  simple  if  the  federation  variation  in  the  structure  has  been  determined. 

For  practical  reasons,  it  is  desirable  and.  Indeed,  often  necessary  to  idealize  the  pro- 
blem to  an  extent  that  ihere  are  minor  effects  superimposed  upon  a major  effect.  These 
minor  effects  are  then  ignored. 

The  stress  measurements  in  these  experiments  were  taken  at  the  mid-snan 
section  so  as  to  minimize  the  effects  of  the  root  and  the  free  end  of  the  cantilever  beam. 
The  idealized  analytical  structure  is  an  infinitely  long  beam  which,  of  course,  cannot 
Iw  attained!  in  the  laboratory.  Instead.  St.  Venant's  Principle  must  be  invoked. 

The  gaaes  were  oriented  so  as  to  measure  stresses  in  the  spanwise  direction 
and,  hence,  chord  wise  stresses  were  i§;iored.  in  addition,  the  gages  were  placed  on 
only  one  side  of  the  skin.  This  procedure  assumes  that  there  is  no  stress  (prodient  across 
the  skin  thickness.  The  assumptions  which  have  been  described  here  or  * those  made  for 
the  elementary  theory  of  the  bending  of  beams. 

The  analysis  of  the  experimental  results  will  be  considered  from  the  follow- 
ing three  viewpoints: 


4. 1 Elementary  beam  theory 

4.2  Distortion  of  the  cress  section 
4.  3 Spanwise  effects. 


4. 1 Elementary  Beam  Theory 

Timoshenko  and  Gccdier  (Ref.  18)  describe  the  physical  reasoning  involved 
in  the  deduction  of  a simple  theory  about  thermo!  stresses.  In  this  simplified  theory  as 
presented  for  application  to  the  experimental  results.  It  is  necessary  to  assume:  (a)  that 
the  temperature  distribution  is  symmetrical  about  the  Z-axIs  (see  Figure  4. 1)  and  (b) 
that  the  temperature  distribution  is  not  a function  of  the  spartwise  coordinate.  Thus, 
the  change  In  temperature  can  br  written  os 

AT  » AT(Sj 

where  5 is  a peripheral  coordinate. 
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The  behavior  of  the  beam  In  adjusting  ihi  distortion  and  internal  stresses  to  a certain 
temperature  distribution  is  explained  as  follows: 

(1)  Each  element  of  the  beam  is  allowed  to  expand  fireel/  as  the  cross 
section  attains  some  particular  temperature  distribution  above  the  ambient  value.  At 
this  stage,  the  beam  is  regarded  essentially  as  an  assemblage  of  independent  elements 
which  were  initially  fitted  together  in  the  shape  of  a beam. 

(2)  In  order  to  fit  the  elements  together  again,  stresses  sufficient  to  com- 
press the  elements  back  to  their  original  sizes  and  shapes  are  applied.  As  In  the  ele- 
mentary beam  theory,  it  is  assumed  that  the  important  distortions  and  stresses  ore  those 
in  the  spanwise  direction.  Hence,  only  the  spanwise  temperature  distortions  are  con- 
sidered, and  the  important  compressive  stresses  which  need  to  be  applied  are  those  in 
the  spanwise  direction.  These  stresses  are  denoted  as  follows: 


tv  fc  \ 

'-A,  i w 


rM-5  ati'c! 

' / W 4—4  - 


where 


E Is  Young's  modulus  which  is  a function  of  the  temperature, 

O?  is  the  coefficient  of  thermal  expansion  which  is  a function  of  the 
temperature, 

AT  Is  the  temperature  rise  above  the  ambient  value. 

(3)  The  beam  con  now  be  re -assembled.  hut  the  force  equilibrium  of  the 
beam  has  been  disturbed  by  the  application  of  the  compressive  stresses  described  above. 
To  restore  equilibrium,  a fictitious  axial  load  and  bending  moment,  statically  eau'va- 
ient  but  opposite  in  character  to  the  summation  and  moment  of  the  compressive  stresses 
across  the  cross  section,  must  be  applied  to  the  free  end  of  the  beam.  On  the  basis  of 
St.  Venant's  Principle,  the  local  effects  of  the  force  and  moment  are  deemed  negligible 
at  cross  sections  not  in  the  immediate  vicinity  of  the  points  of  application. 

(4)  The  stress  distribution  at  cross  sections  far  removed  from  the  roct  ardl 
the  free  end  is  obtained  by  tho  superposition  of  the  effects  described  in  (2)  and  (3)  above 
as  follows: 


where 


(4.2) 

(4.3) 


fsl  '*^1,  'T*  V 


fi 


<y  - - 

is 

F is  the  statically  equivalent  axial  load  of  opposite  character 


(4.4) 


M is  the  statically  equivalent  bending  moment  of  apposite  character 
A is  the  cross- soctiooal  area 

I«l  is  the  moment  of  inertia  of  cross  section  about  the  y-axls. 

The  axial  load,  F , and  bend’  i»y  mOmuTii;  M , applied  at  the  free 
end,  are  obtained  from  the  tallowing  expressions: 


F^-Jc^dA  = j EaATdA 

A A 

M -[-fo  ZAA]  • - 1 E 5 ATZ  JA 


(4.5) 


(4.6) 


For  purposes  of  calculation,  it  is  necessary  to  replace  the  integral  over 
the  cross  section  by  a summation  because  the  experimentally  determined  temperature 
distribution,  AT (S),  cannot  be  written  easily  in  analytical  form.  The  cross  section  of 
the  experimental  box  beam  was  divided  into  forty-two  elements  (see  Figure  4. 2)=  Only 
twenty-one  baa  to  be  considered  because  of  the  model  symmetry  and  the  assumed  tem- 
perature symmetry  about  the  2- axis.  The  expressions  for  the  three  components  of  the 
total  stress  are 


(q;<1)i=-E«ATi 

<^A‘ii2EXATi 


(4.7) 


(4.9) 


r 2 2’e**ta  (4-9) 

“3  i«l 

vfhere  the  subscript,  l , denotes  each  element.  The  values  for  of  and  E were  corrected 
?«r  temperature  a*  shown  in  Figure  4. 3. 

The  stresses  as  calculated  from  the  above  expressions  are  tabulated  in  Table 
A.  4..  It  must  be  remembered  that  the  temperatures  used  in  these  calculations  were  ex- 
trapolated and  interpolated  from,  eleven  meej^rad  tesperotyres: 
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These  calculated  stresses  together  with  the  experimental  values  are  shown 
for  heating  times  of  30,  60,  90  and  120  seconds  ?n  Figures  4.  4,  4. 5,  4. 6 and  4. 7 for 
Configuration  1 and  in  Figures  4.3,  4.9,  4.  10  and  4. 11  for  Configuration  2.  The 
agreement  in  the  bottom  skin  is  fair  while  the  agreement  in  the  web  is  poor. 


4. 2 The  Distortion  of  the  Mid-Span  Cross  Section 

If  the  box  beam  is  visualized  as  a collection  of  closed  rigid  frames  (see 
Figure  4. 12),  then  the  existence  of  a temperature  change  can  be  seen  to  cause  a dis- 
tortion of  the  rigid  frame.  Since  the  rigid  frame  is  an  indeterminate  structure,  this 
distortion  v/III  induce  both  axial  stresses  and  bending  stresses,  in  this  case,  the  terms, 
"axial  stresses"  and  "bending  stresses",  are  to  be  considered  with  reference  to  the  be- 
havior of  the  rigid  frame. 

The  concept  of  an  elastic  center  together  with  the  superposition  equations 
will  be  used  in  this  analysis  (see  Ref.  14).  These  superposition  equations  ore  applied 
to  the  primary  structure  (see  Fig.  4. 13)  which  is  statically  determinate. 


daa 

$a.b 

^ba 

Sbb 

^ca 

Scb 

where 


= 0 <*U>) 


Xa „ Xc,  are  the  redundant  bending  moment,  axial  load  and  shear 

force  which  act  at  the  elastic  center  (see  Figure  4. 13) 

6 mn  is  the  relative  displacement  of  the  points  of  application 
of  redundant  m due  to  a unit  value  of  redundant  n 

S rti  T is  the  relative  displacement  of  the  points  of  application 

of  redundant  TO  due  to  the  temperature  distribution,AT(^. 
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FIGURE  4. 13  PRIMARY  STRUCTURE 


It  will  be  assumed  that  the  temperature  Is  uniform  across  the  thickness  of 
the  frame.  The  virtual  work  expressions  for  the  Influence  coefficients.  6 and 
6mT,  are  as  follows: 
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S»o  "2  jj  M*  M<j  ex'  (b  £AJ 

♦JV£ } 

Sab  -iff a Mb  ff  +/k  G. 

Sic  -ifKM.If  *fFa  £ 

Sir=2«Jfi  ATds 


SbT=5acjFb  zrrds 


Scr^Jfi  AT ds 


whore 


E Is  Young's  modulus 
▼ • .i-  - i .t  i i:„  _ ^jjj: 

X IMC  I9IUIII9III  V*  ii*9t  ••%■  — 

A Is  the  cross-sectional  area  « (l)t 

Mm  is  the  bending  moment  in  the  primary  structure  caused  by  a unit 
value  of  redundant  vn 

Fwt  1$  the  axial  load  in  the  primary  structure  caused  by  a unit  value  of 
redundant  m 

Mo  Is  the  bending  moment  in  the  primary  structure  caused  by  external 
loads 

Fo  is  the  axial  load  In  the  primary  structure  caused  b>  external  loads. 
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From  the  symmetry  of  the  cross  section  and  of  tho  temperature  distribution, 
the  following  results  are  obtained: 
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(4.13) 


The  temperature  distribution  of  Table  A.  2 for  a heating  time  of  120  seconds 
has  been  used  to  evaluate  Sbj  numerically  with  the  result  that 


E Sb7  - 5^ Fb^S  AT- 529.6 oc  * -6880 *IO~s* 
E Sbb  - 5. 93  x io"4  + 0. 1 32  * IO"4-  6.06 x (O'4 
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Tli©  redundant  axial  load,.  Xj,  , 


has  the  value 


y,  =-~6»,3g3*IQl-  g i|  315  ik5 

Ab  6l06 * iO~*  11.07  IDS. 


{4. 15) 


ure  4. 14. 


Plots  of  the  bending  moment  and  axial  load  distribution  are  shawm  In  Flg- 


17  in  ee  I 


-17  in# 


+ I7in# 


17  in# 

V 


“17  in* 


FIGURE  4. 14  BENDING  MOMENT  DISTRIBUTION 


At  point  <S  in  Figure  4. 14,  the  maximum  stresses  are 


(a-)  =_£ji±-E. 

lus/e  I A 


The  e*jl*ide  fiber  at  point  € has  the  stress 
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(CTs/e  (|)(|2XQ274)(!0r5  + (0(0.065; 


(4.17) 


= 7 


pst 


The  inside  fiber  at  point  e has  the  stress 

(Oi)e  =-l68-l75  = -343Psi  M-'S) 

Thus,  the  stresses  which  result  from  the  distortion  of  the  cross  section  are  inconse- 

»e_| 

quomiuis 


4, 3 Sponyylse  Pffach 


The  principle  of  minimum  potential  energy  can  be  used  to  obtain  approxi- 
mate solutions  to  the  three-dimensional  equations  of  elasticity.  If  there  are  no  external 
or  body  forces  acting  upon  the  structure,  the  differential  equations  of  equilibrium  have 
the  following  form  (Ref.  18): 

(•>'■'<5)  ax  (ft  §if  )■*■<=  V'U-  pfv  Tlf-'O 

^*s)  % (ft-  V v - Wv  Is  * o <“•  ”) 

(A  *e)  ft  (ft If ) *<s  v1  w - Jg,  H • o 


where  the  constants  are  defined  as 

F 


X = 


(i+j/’Xi-zy) 


2 ( i + ir) 
if  - Poisson's  ratio. 
The  stress-strain  relations  become 
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When  properly  minimized,,  the  potential  energy  expression  which  w!!E  yield 
the  above  equations  of  equilibrium  and  stress-strain  relations  has  the  form  shown  below: 

kJJjf  + + €z)  +2-G fee*4 £'j  4 €l  ) 

(4. 22) 

+6(«JS  V V€*)}d*  J*  dz 


22) 


It  is  readily  shown  that  the  vanishing  of  the  variation  of  ths  above  integral  will  lead  to 
the  three  differential  equations  of  equilibrium.  Furthermore,  If  the  strain  energy  density 
function  Is  defined  as 


W {A(e* 4 e3  4 €if+  2s(e/+ €/♦€/)  *G  fa,* fa,** JJ,*) 


(4.23) 


!-2i7 

it  is  also  readily  shown  that 
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(4.24) 


ASRL  TR  25-13 


-'^"ScV1-3 


**  '■'***’  ^ 


Sf  4.  . 


Alfl&SsasHsa 


Thus,  It  Is  scan  that  the  potential  energy  function  satisfies  all  the  prerequisites  neces- 
sary for  the  ihfss-s:n»rj!cno!  e^uchcns  of  thermos  lastlclty  in  which  external  and 
body  forces  are  absent. 

4-.  3.  V Application  to  the  Box  Beam 

The  principle  of  minimum  potential  energy  as  stated  will  be  specialized  to 
the  problem  of  the  box  beam.  It  will  be  necessary  to  make  the  following  assumptions: 

(1)  The  cross-sectional  shape  is  preserved  by  many  closed  diaphragms  which 
are  rigid  in  their  own  plane  but  completely  flexible  out  of  their  plane,  he.,  there  is 
no  resistance  to  warping.  In  addition,  it  is  necessary  to  assume  that  these  diaphragms 
are  perfectly  insulated  from  the  rest  cf  the  structure. 

(2)  The  wall  (or  shell)  or  the  box  beam  is  in  a state  of  membrane  stress, 

I.e.,  the  stresses  and  hence  the  temperature  is  com  ram  aver  in*  mibRness  cf  the  she!!. 

These  assumptions  are  the  far-: liar  ones  used  in  the  analysis  of  the  effects  of  warping 
restraint  on  the  stress  distribution  in  sneii  beum»  (Refs.  !!,  !2  and  13)- 

The  expression  for  the  potential  energy  of  strain  is  thus  reduced  to  the  fol- 
lowing expression: 
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t 

1 

I' 

i 


-2E' "dc -AT  {£,+€*)]  JsJx  (4.25) 


The  pertinent  stress-strain  relations  are 

C5^  = E ' (e+  E *oc  n~T 


where 


<ts  - E'(es-  ye,)- 

E'‘  ;rf* 

e'.-E. 

t I- V 

r _ _ E 
6 ' *0*-.-) 

)'  is  the  shear  strain 

£x  is  the  strain  in  a spanwise  direction 
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€s  is  the  strain  In  a peripheral  direction 

is  the  normal  stress  in  a spanwise  direction 
C%  Is  the  normal  stress  In  a peripheral  direction 
of  is  the  coefficient  of  thermal  expansion 
AT  is  the  change  in  temperature 
t Is  the  thickness. 

Points  on  the  crass  section  are  located  by  a peripheral  coordinate,  s , 
and  an  angle,  oc  , which  the  tangent  ot  point  & makes  with  the  y -axis  (see 
Figure  1. 16).  Displacement*  In  the  ««««.'!*«  f x ).  vertical  ( B ) and  peripheral 
( S ) directions  will  be  denoted  by  u , u/-  and  ^ , respectively.  Rotation- 

al end  horizontal  displacements  are  neglected  since  both  the  cross  section  of  the  model 
and  the  temperature  distribution  are  symmetrica!  about  the  Z -axis.  A vertical  dis- 
placement of  the  cross  section  causes  a peripheral  displacement  of  the  amwni  (m«  Fig- 
ure 15) 

Jf=U/-sinX  (4.28) 


Z 


FIGURE  4. 15  RELATION  BETWEEN  VERTICAL  AND 
PERIPHERAL  DISPLACEMENTS 
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Thus,  the  potential  energy  expression  becomes 

7fA  * 2. ffi  ^E'  Of  * +■  u>tj.  sin  or) -2^Eu*AT*Uxj  dsdx 

where 


~OU 

r>A 


U, 


/* 


■QU  "B  5 

■OS  -ox 


= Us-t-u£ 


sin  OC 


(4.30) 


The  peripheral  strain,  C3  , vanishes  because  the  assumption  of  rigid  diaphragms  pre- 
vents the  occurrence  of  such  strains. 


If  has  been  found  that  the  temperature  distribution  at  any  particular  time 
can  be  adequately  represented  by  an  expression  of  the  form, 

^T(S)=<Ae'^'S(eas^-s+si„5:S/)*B  (4.31) 


where  ~jr  _ 

and  D are  constants 

S0  “ 3/2  total  peripheral  length  = 9 inches. 

Figures  4.  17  and  4. 18  are  plots  which  show  the  agreement  between  the  experimental 
temperature  distribution  and  the  distribution  given  by  the  above  assumed  function. 

This  agreement  suggests  the  following  assumed  function  for  the  spanwise  displacement: 

\J  (x,.s)  (x)  I ( 2S-  9)  sir,  CK-  S cos  « ( +■ 

_ t,  ^ (4-  22) 

jpz  (x)e  cos -s + (y) e 5o  s\„  — - s 

The  first  term  is  not  a function  of  S and,  hence,  represents  a rigid  translational 
displacement  in  which  all  points  of  the  cross  section  move  ihe  same  amount.  The  coef- 
ficient of  $6,  (x)  represents  a rigid  planar  displacement  which  corresponds  to  the 
application  of  a bending  moment  about  the  if  -axis.  The  warping  modes  are  [repre- 
sented by  <^2.  ant^ 

With  the  substitution  of  these  expressions  for  AT  and  U , the  pot- 
ential energy  expression  be  vCm»6S 
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$1 
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TTS=2.‘X  f J t |E  [Vox "*■ 9*tx  (2S  s'mOC-9s»nC(-3 CO SGf) 

"o  “o  *-  L _ 

+>4.^  ‘cosas  + ^x  e sin  asj 
+g[V»  ^ (2S  sin  <X-9 sinoc-3coscx)~ae^S(cGS  as  -*  Sin  as)^ 
+ ae~tcS(cos  as- sin  as)^-*-U/y  sin  <*] 

-2^\e  a5(cos  as  ♦ sin  as)  + B][$4*+^xfe>  si*,cx’~  9 sin cX"  -3cascx) 


'2. 

(4.33) 


d -_as 


00s  as 


+ e 4S  sin  as  |j  ds  d>c 


Wooi‘5 


Jl  U the  span  length 

„ tt 
a-  9 

A E“X 

B * of  E"  B 

It  is  convenient  to  make  the  following  substitutions: 
-p,  (s)  = e'^cos  as 
e-*'*  sin  as 

.p>,(s}«  e"^*J(cas  as  -t  Sin  as)16^ 
_p4(s)=  e'^^c-tts  as  - sin  as)=--p,—-px 
(25-9)  sinoc  - 3 cos  or 

f*(s)= 

svnOC 


(4. 34) 


(4. 35) 
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The  strew -*"ro In  expression  now  becomes 

q;  * £ ' {&«  * f6  f,  * ] - [a?**  b]  <«•  34) 


The  variation  of  the  Integral  yields 

S1Ts-2^  Sfaf.SfeiiSfa] 


+ S^xj|  ds  d/ 

j 

It  Is  noted  that  the  following  integrals 

appear  In  the  above  expression: 
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Clj  « zfS* t^*d3-  0/651 
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art’  2^t^*ds  = 15.42/ 
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(X.yt~z£  ty»xf,^s=-.,.l/34  (4.38) 
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a5  » zj^' t^f  ^ ds  = 0.07455 

a8 =-0.5498 
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These  constants  were  evaluated  on  the  basis  of  the  Idealization  of  the  cross  section 
shown  in  Figure  4. 19. 

The  vanishing  of  the  variation  of  the  strain  energy  yields  the  following 
differential  equations: 


9^oxx  ^ 

- e ' (a,9  ^ +al5  fa*  ■ o w 

- E (a, ft**  +a* p,**  9w+<Z,  (c)  "3?> 

“ t (C<2  JPqXX  PlM  P'zxx  "’Ui4  /^AX  / " “V  ) C‘— 

+ (e) 


The  natural  boundary  co'td’ti  ~*i$  ar  A-Jt  ore 


(c)  (4.40) 


E'fa0j<„  *a,^u*a,Ax)  -a  (a,  ♦a*)*  b<2„ 

E ($13  ^IX  +^hS^2J<  +^!5^fx/  =A(l2H  +<^|5/ 
E'('a,$w+<2rt  ^.«*a3  A*  *a,*u)‘A  (a3+as)+B0., 
t ’ {Oi^c-t  rG ;g  ^ix  -a,  j4x  “ Aids+O.j)*'  B&, 

-o 


From  differential  equation  (4. 39e)  and  boundary  condition  (4.  4Ge),  the 
fo'lowhtg  is  true  everywhere 
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FIGURE  4. 19  IDEALIZATION  OF  CROSS  SECTION  FOR  CALCULATION  OF  <Z,  , 
az  , . . . al9  PARAMETERS 
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(4.41) 


a, 


^■i» 

It  can  also  be  shown  that 


alo*  N 


a,fe  -a  !9  ” (<** - * 


14. 42) 


and  that 


(4.43) 


(4.44) 


• hi*  rosu'ir  reduces  equation  oro;  to  me  :c: lowing  rorm: 

— E f<3.  J3  14  +^jsr  ^<5 xx V = v^‘ 

Thus,  equations  (4. 39a)  and  (4. 39b)  together  with  boundary  condition:  . 
(4. 40a)  and  (4. 4Cb)  insure  the  satisfaction  of  axial  force  and  bending  moment  equili- 
brium at  every  cross  section.  These  relations  are 

e' (Oo -k2z - [a  (a  + Ba0  ] =o  (4.45) 


E (&»3  fizx  +G-I9  5^4x)  " [a  (^14  )]  sO 


(4.46) 


By  eliminating  <floxx  and  from  the  equctions  (4. 45)  and  (4. 46),  the  two  dif- 

ferential equations  for  the  warping  modes,  %bI}  and  ^ , are  seen  to  be 


-E  ’ fcaj-  ) &x*  * (a*  " 

.g{ 


a,ai  a.«a,r 


a, 5 / /T**X 
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>6  {-«;&■*  <2b^)‘0 


(4.43) 
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The  solutions  for  £z£r  and  ^ are  of  the  form 

*=■  Fj.^  sinh  rjX  + sink  l^x 

*£  * £(0  smh  r,X  ♦ R W -Slnh  rxX 


(4.49) 


where 

r,  =0.56636 

(4.50) 

t£  = 1.63930 

The  hyperbolic  sine  terms  must  vanish  in  order  io  satisfy  the  boundary  conditions  at  the 
root  of  the  beam. 


ir.  order  to  illustrate  the  spnnw?**  effects,  the  constants  for  the  temperature 
change  in  equation  (4. 31)  were  taken  to  be  (see  Figure  4. 18) 

A - 216 
B =38 

These  eon!*"nts  approximate  the  experimental  temperature  distribution  at  the  end  of 
the  heating  period  of  120  seconds.  The  remaining  constants  for  0m  and  y>^  were 
obtained  from  boundary  conditions  (4. 40c)  and  (4. 40d).  Thus, 


rr-1 

0.9585*10'*  0.1335*  IO*9 

Lr  J = 

-OA782xtO'e  0.4761*10'* 

(4-  51) 


By  substituting  equations  (4. 49)  and  (4-  51)  into  equation  (4. 36)  me  normal 
stress  expression  becomes 


(4.52) 


St  car.  be  shown  that  if  the  span,  Jl  , is  infinitely  long,  the  warping  modes, 
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^ will  disappear.  Thus, 


-p'  + 


^15 


(4.53) 


A good  measure  of  the  effect  of  the  free  end  is  seen  to  be  the  two  expressions. 


lA-r'  cC,l 

U ' J (4.54) 

[a-eV4x] 

These  are  shown  plotted  in  Figure  4. 2C.  The  effects  of  the  fires  end  are  thus  shown  to 
be  reduced  to  negligible  proportions  at  a distance  of  about  one-third  of  the  total  peri- 
pheral length  of  18  inches  from  the  free  end. 

Figures  4.  21,  4. 22,  4. 23  and  4. 24  are  comparisons  of  the  stress  distribu- 
tions corre^or-ding  to  the  above  analysis  and  the  experimental  stresses  at  heating  times 
of  30,  60,  90  and  120  seconds,  respectively.  The  agreement  is  reasonably  good  for 
30  and  60  seconds  but  rather  poor  ot  90  and  1 20  seconds.  The  latter  discrepancies  may 
be  uue  to  the  thermal  buckling  of  the  lower  skin  which  would  change  the  stress  distri- 
bution. 


4. 4 The  Deflection  at  the  Free  End 


The  deflection  at  the  free  end  car.  be  calculated  accuraTely  bv  6 temennsry 
beam  theory.  From  this  theory,  it  can  be  shown  that  the  tip  deflection  of  a cantilever 
beam  acted  upon  by  a tip  bending  moment  is  given  by 

S-Hf  «.55> 


Thw  analyses  c?  Section  4. 1 yielded  bending  moments  of  -20, 970  in.  -lbs.  and  -20,  360 
in.  -lbs.  for  Configurations  1 ana  2,  respectively,  at  the  end  of  the  heating  period. 
Thu*, 
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FIGURE  4.21.  COMPARISON  OF  THEORETICAL  AND  EXPERIMENTAL  STRESS  FOR 
30  SECONDS  - CONFIGURATION  2 

(Note:  Stresses  are  plotted  normal  to  the  center  line.  Compressive  stresses  are  plotted 
on  the  outside  of  the  cross  section,  tensile  stresses  on  the  inside.) 
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FIGURE  4.22  COMPARISON  OF  THEORETICAL  AND  EXPERIMENTAL  STRESSES 
FOR  60  SECONDS  - CONFIGURATION  2 

(Note:  Stresses  ore  plotted  normal  to  the  canter  line.  Compressive  stresses  ore  plotted 
on  the  outside  of  the  cross  section,  tensile  stresses  on  the  Inside.) 
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FIGURE  4.24  COMPARISON  OF  THEORETICAL  AND  EXPERIMENTAL  STRESSES 
FOR  120  SECONDS  - CONFIGURATION  2 


{Note:  Stresses  are  plotted  normal  to  the  center  line.  Compressive  stresses  are  plotted 
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tensile  stresses  on  the  inside. ) 
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These  calculated  values  compare  favorably  with  the  observed  displacements  of  0. 168 
Inches  and  0. 161  Inchss,  respectively. 
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SECTION  V 


CONCLUSIONS  AND  RECOMMENDATIONS 


(1)  The  primary  mode  of  heat  transfer  In  these  experiments  was  conduction. 
Convection  and  radiation  were  relatively  unimportant  because  of  the  precautions  taken. 

(2)  The  agreement  between  the  experimental  stress  distribution  end  the  stress 

distributions  calculated  by  two  different  methods  on  the  basis  of  temparstures  meai-.T«d 
at  eleven  points  (see  Sections  4. 1 and  4.  3)  is  fair  for  heating  times  of  30  and  60  seconds. 
But  the  agreement  is  poor  for  times  of  90  and  120  seconds.  This  variation  suggests  that 
buckling  of  the  lower  skin  may  have  occurred  although  there  was  no  discernible  evidence 
at  the  conclusion  of  tho  rests,  it  was  «*;«  poi'/.bie  iz  cbsrrve  the  surface  during 

the  tests,  and  the  buckling  of  the  skin  could  account  for  the  lower  experimental  stresses 
at  these  two  longer  heating  times.  Another  reason  for  the  differences  at  all  times  is  that 
the  cunuiiioTu  ot  the  corner  fhrngei  co*«ld  n->t  be  determined  experimentally  with  sufficient 
detail  as  to  temperatures  and  stresses.  The  relatively  large  flange  areas  play  a dominant 
role  in  the  satisfaction  of  force  and  moment  equilibria  c«er  the  cross  section.  Thus,  an 
error  in  the  determination  of  the  flange  temperature  will  be  magnified  In  the  calculations 
of  the  stresses  aver  the  entire  cross  section. 

(3}  The  highest  stresses  were  in  the  region  of  the  highest  temperatures.  Thus, 
the  largest  stresses  were  In  the  bottom  skin  which  received  the  direct  thermal  radiation 
and  were  compressive  In  nature.  The  smallest  stresses  occurred  lr»  the  upper  skin. 

(4)  As  expected,  the  center  portion  of  the  vertical  web  was  in  tension  In 
order  to  counterbalance  the  compressive  stresses  In  the  bottom  skin. 

(5)  By  contrast,  the  tension  in  the  lower  legs  of  the  corner  flanges  on  the 
bottom  of  the  cross  section  was  not  expected.  Measurements  at  the  bottom  of  the  corner 
flange  could  not  be  obtained  because  of  the  difficulty  of  squeezing  a strain  gage  between 
bolt  holes.  The  calculated  stresses,  however,  indicate  tension.  Thus,  the  behavior  of 
the  beam  is  such  that  the  tensile  stresses  In  the  flanges  on  the  bottom  help  to  balance  out 
the  bending  moment  Induced  by  the  large  conpressive  stresses  In  the  heated  skin.  Inis 
effect  relieves  the  load  which  the  upper  skin  would  have  hod  to  carry  if  the  entire  bot- 
tom of  the  beam  had  been  in  compression  and  explains  why  such  low  stre»$  levels  were 
encountered  in  the  unheated  upper  surfocn. 

(6)  The  severest  temperature  gradient  occurred  In  the  bottom  and,  hence, 
was  In  the  chordwise  direction.  It  is  suspected  that  the  furnace  radiates  ci  flux  which  is 
not  uniformly  distributed  over  the  cross  section  end  which  foils  off  rather  significantly  at 
the  edges.  Th?«  uneven  distribution  contributes  to  a chordwise  temperature  gradient.  But 
even  with  a uniformly  distributed  therms!  flux,  the  flow  of  heat  into  the  relatively  large 

ASRL  TR  25-1 3 - 62  - 


masses  represented  by  the  corner  flange:  and  from  them  Into  the  web  would  cause  a 
rather  large  chordwlse  gradient. 

(7)  The  distortion  of  the  cross  section  due  to  the  Increase  In  'emperature 
Induces  negligible  stresses  In  the  chordwlse  direction. 

(8)  The  effect  of  the  free  end  on  the  stress  distribution  Is  confined  to  a 
region  within  one-tnlrd  of  the  total  peripheral  length  from  the  tip.  Since  the  experi- 
mental stresses  were  measured  at  mid-span  which  is  equal  to  two-thirds  of  the  peripheral 
length,  they  were  unaffected  by  ths  finite  span.  Thus,  the  model  was  long  enough 

to  bs  considered  on  ?of?nltu  beam. 

(9)  It  Is  felt  that  the  conditions  at  the  corners  should  be  investigated  In 
greater  detail  despite  the  difficult  experimental  problems. 

(10)  The  characteristics  of  the  radiant  furnace  should  also  be  determined. 

A calculation  of  the  temperatures  in  the  structure  presumes  a knowledge  of  the  heat 
source,  and  the  electrical  power  Input  to  the  heater  Is  not  sufficient  data  for  such  com- 
putations. 


(1 1)  The  high  compressive  stress  levels  In  the  skin  as  compered  to  the  low 
streu  levels  In  the  flanges  suggest  that  thermal  buckling  will  be  a serious  problem. 
These  large  skin  stresses  occur  at  fairly  high  temperatures  so  that  thermoplastic  rather 
than  thermoelastic  buckling  will  need  to  be  Investigated. 
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PHYSICAL  PROPERTIES  OF  24  S-T  AND  61  S-T  ALUMINUM  ALLOY 


PROPERTY 

24  S-T 

61  S-T 

g 

Specific  Gravity,  gny^cm 

2.77 

2.70 

Weight,  lbs/in^ 

0.100 

0. 098 

Ultimate  Tensile  Stress,  lbs/in^,  at  60°F 

68000 

45000 

2 1 2“? 

: 1 AAA 

0 1 VW 

-r  i ww 

300°F 

431/00 

32000 

400°  F 

26000 

i?C0C 

500°F 

14000 

7000 

Yield  Tensile  Stress,  lbs/in^  at  68°F 

48000 

40000 

212°F 

45000 

37000 

nonop 

37000 

30000 

400°F 

22000 

16000 

500°F 

10000 

5000 

E longation  in  2 in,  % 

19 

12 

A 

Modulus  of  Elasticity,  Ibs/in* 

io.5  x io6 

10.1  x 10^ 

Approximate  Melting  Range,  °F 

935-1180 

1080-1205 

Plastic  Flow  Range,  °F 

600-900 

600-900 

Specific  Heat,  BTU/!u/°F 

0. 23 

0.22 

Coefficient  of  Thermal  Conductivity, 

BTU/hr/ft2/°F/ft,  (0-212°F) 

70.  1 

89.5 

Average  Coefficient  of  Thermal  Expansion 

/°F,  68-2 12°F 

12.9  x 10"6 

13.1  x 10~° 

6 8 -392° F 

13.3  x 10"* 

13.5  x 10"® 

68-572°F 

13.7  x 10"® 

14.1  x lO^6 

Modulus  of  elasticity  arsd  coefficient  of  thermal  expansion  vs.  temperaruie  are  shown 
in  Figure  4.  3 and  were  obtained  from  References  1 and  2. 

The  above  physical  properties  have  been  compiled  from  References  1,  2 , 3,  11,  and  12. 
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TABLE  A.  2 


EXPERIMENTAL  TEMPERATURE  VARIATION  AT  EACH  GAGE  LOCATION/  °F 


GAGE 

NO.  T-l 

G AGE  NO.  T-2 

GAGE  NO. 

T-3 

CONFIGURATION 

CONFIGUR 

ATiON 

CONFIGURATION 

TIME 

SEC. 

1 

2 

1 

2 

1 

2 

0 

78 

78 

/O 

TO 

• W 

78 

78 

10 

78 

78 

78 

78 

78 

78 

20 

78 

78 

78 

78 

79 

78 

30 

79 

78 

79 

79 

80 

80 

40 

81 

80 

81 

wu 

82 

81 

50 

82 

80 

83 

81 

86 

83 

60 

85 

82 

86 

83 

88 

85 

70 

88 

84 

90 

95 

93 

89 

80 

92 

86 

94 

89 

99 

93 

90 

97 

91 

100 

92 

104 

97 

ICO 

103 

95 

105 

97 

111 

103 

110 

110 

100 

no 

101 

117 

108 

120 

115 

104 

117 

106 

124 

114 

m 


TAB!?  A-  2 (Continued) 


EXPERIMENTAL  TEMPERATURE  VARIATION  AT  EACH  GAGE  LOCATION,  °F 


GAGE  NO.  T-4 

GAGE  NO. 

T-5 

GAGE  NO. 

T**6 

CONFIGURATION 

CONFIGURATION 

CONFIGURATION 

TIME 

SEC. 

1 

2 

1 

A. 

1 

2 

0 

78 

78 

■70 

• o 

78 

78 

78 

10 

79 

79 

81 

80 

R3 

81 

20 

M "• 

w • 

Ml 

BU 

84 

88 

86 

30 

84 

84 

89 

88 

95 

92 

40 

88 

87 

95 

93 

103 

99 

50 

93 

91 

102 

99 

112 

106 

60 

99 

96 

109 

105 

121 

114 

70 

105 

100 

116 

111 

130 

122 

80 

• • 
i \t 

1 i\f 

1 W 

124 

na 

140 

131 

rV\ 

vu 

119 

112 

133 

125 

150 

140 

10C 

127 

119 

142 

133 

161 

148 

no 

135 

125 

151 

141 

170 

157 

120 

144 

133 

160 

149 

180 

167 
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TABLE  A.  2 {Cor*tniM»d) 


EXPERIMENTAL  TEMPERATURE  VARIATION  AT  EACH  GAGE  LOCATION-  °F 


GAGE  NO.  T-7 

GAGE  NO. 

T-8 

GAGE  NO. 

T-9 

CONFIGURATION 

CONFIGURATION 

CONFIGURATION 

TIME 

SEC. 

1 

2 

1 

2 

1 

1 

9 

6m 

0 

78 

78 

78 

78 

78 

78 

10 

84 

82 

86 

83 

108 

103 

2V 

93 

7U 

97 

91 

129 

123 

30 

103 

99 

108 

101 

148 

141 

40 

113 

109 

120 

112 

166 

159 

50 

123 

118 

132 

122 

183 

174 

60 

135 

127 

145 

134 

199 

190 

70 

146 

140 

157 

145 

124 

205 

80 

157 

150 

170 

156 

229 

218 

90 

168 

160 

182 

167 

243 

233 

100 

179 

170 

195 

179 

259 

247 

110 

189 

181 

206 

189 

272 

259 

• AV 

199 

191 

218 

201 

285 

272 
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TABLE  A.  2 (Ccntlnm*!) 


EXPERIMENTAL  TEMPERATURE  VARIATION  AT  EACH  GAGE  LOCATION,  °F 


TIME 

GAGE  NO.  T-10 
CONFIGURATION 

GAGE  NO.  T~!l 
CONFIGURATION 

* 

i 

2 

1 

2 

A 

V 

78 

73 

78 

78 

10 

116 

110 

120 

111 

20 

145 

140 

149 

143 

30 

169 

163 

177 

171 

40 

191 

185 

202 

194 

50 

210 

203 

222 

216 

60 

229 

223 

242 

235 

70 

246 

240 

260 

254 

80 

262 

255 

279 

271 

90 

277 

269 

294 

288 

100 

294 

287 

312 

3C3 

110 

307 

300 

326 

318 

120 

322 

314 

341 

332 
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▼ ABIC  A T 


9 HVM  v* 

EXPERIMENTAL  STRAIN  AND  STRESS 


STRAIN  GAGE  NUMBER  S - 1 


TIME 

CONFIGURATION 

1 

CONFIGURATION  2 

SEC-. 

«•  10 

6 E- 10-6 

(P*0 

e-106 

E* 10“^ 

(psl) 

10 

- 2 

10.5 

- 21 

- 4 

10.5 

■ 42 

20 

8 

10.5 

84 

1 

10. 5 

i 1 
1 1 

30 

15 

10.5 

160 

15 

10.5 

120 

40 

26 

10.5 

270 

25 

10.5 

260 

50 

33 

10.5 

350 

31 

10.5 

330 

60 

27 

10.5 

280 

51 

10.5 

540 

70 

19 

10.5 

200 

64 

10.5 

670 

80 

18 

10.49 

190 

53 

10.5 

560 

90 

8 

10.49 

80 

44 

10.48 

460 

100 

- 5 

10.48 

- 50 

41 

10.47 

430 

110 

-15 

10.47 

-160 

29 

10.46 

300 

120 

-15 

10.46 

-160 

29 

10.46 

300 
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EXPERIMENTAL  STRAIN  AND  STRESS 
STRAIN  GAGE  NUMBER  S - 3 


TIME 

CONFIGURATION  1 

CONFIGURATION  2 

SEC. 

e.  106 

E-10"6 

(psl) 

.«6 
«•  IU 

E10"6 

(psl) 

10 

- 7 

10.5 

- 74 

- 8 

10.5 

- 84 

20 

- 7 

10.6 

- 74 

- O 

i r\  e 
IV.  V 

U A 

30 

- 10 

10.5 

- 105 

- 20 

10.5 

- 210 

40 

- 27 

10.5 

- 280 

- 37 

10.5 

- 390 

50 

- 33 

10.5 

- 350 

- 54 

10.5 

- 570 

60 

- 64 

10.5 

- 670 

- 65 

10  5 

- 680 

70 

-122 

10.49 

-1280 

-105 

10.5 

-1100 

80 

-130 

10.49 

-1360 

-114 

10.49 

-1200 

90 

-150 

10.48 

-1540 

-120 

10.48 

-1260 

VOO 

-168 

10.47 

-1760 

-155 

10.47 

-1620 

110 

-182 

10.46 

-1900 

-176 

10.46 

-1840 

120 

-203 

10.45 

-2120 

-187 

10.46 

-1960 

I 

I 

I 
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TABlic  A.  3 (Continued) 
EXPERIMENTAL  STRAIN  AND  STRESS 


STRAIN  GAGE  NUMBER  S - 4 


TIME 

CONFIGURATION  1 

CONFIGURATION  2 

SEC. 

e-106 

E-10-6 

(psl) 

e-106 

E-10-6 

(psi) 

10 

23 

10.5 

245 

14 

10.5 

140 

oo 

vv 

in  s 

350 

24 

10,5 

250 

30 

40 

10.5 

420 

31 

10  5 

330 

40 

47 

10.49 

490 

Or 

in  40 

■ »*-  • • - 

410 

cn 

<J>V 

52 

10.49 

540 

43 

10.49 

450 

60 

53 

10.48 

560 

46 

10.48 

480 

70 

57 

10.47 

590 

49 

10.47 

520 

80 

58 

10.46 

610 

51 

10.46 

540 

90 

60 

10.46 

630 

51 

10.46 

540 

100 

60 

10.45 

630 

50 

1C.  45 

520 

no 

57 

10.44 

590 

49 

10.45 

5l0 

120 

57 

10.43 

590 

47 

10.45 

490 

-74- 


ASRL  TV  25-13 


TABLE  A.  3 (Continued) 


EXPERIMENTAL  STRAIN  AND  STRESS 
STRAIN  GAGE  NUMBER  S - 5 


TIME  CONFIGURATION  1 CONFIGURATION  2 


SEC. 

e-106 

5- 10"6 

(psi) 

e.106 

E-10"6 

(psl) 

10 

57 

10.5 

600 

46 

10.5 

490 

20 

104 

10.5 

1090 

86 

10.5 

9«G 

30 

140 

10.49 

1470 

115 

10.5 

1210 

40 

164 

10.48 

1720 

14! 

10.49 

1480 

50 

156 

10.48 

1950 

160 

10.48 

1680 

60 

206 

10.47 

2160 

174 

10.47 

1820 

70 

218 

10.46 

2280 

189 

10.47 

1980 

S3 

227 

10.45 

2370 

200 

10.47 

2100 

90 

234 

10.44 

2440 

204 

10.45 

2130 

100 

257 

10.43 

2680 

211 

10.45 

2210 

no 

265 

10.42 

2760 

214 

10.43 

2230 

120 

268 

10.40 

2790 

221 

10.42 

2310 
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TABLE  A.  3 (Continued) 


EXPERIMENTAL  STRAIN  AND  STRESS 


STRAIN  GAGE  NUMBER  S - 6 

TIME  CONFIGURATION  1 CONFIGURATION  2 


SEC. 

e-  I06 

E-10"6 

(P*0 

e*  106 

E-10"6 

(psl) 

10 

73 

10.5 

760 

-»p 

/J 

10.5 

790 

20 

115 

10.49 

1210 

120 

10.5 

1260 

30 

151 

10.48 

1590 

164 

10,49 

161 0 

40 

176 

10.47 

1850 

131 

10.48 

1900 

50 

193 

10.46 

2010 

199 

10.47 

2080 

60 

206 

10.45 

2160 

210 

10.46 

2200 

70 

214 

JO.  45 

2230 

216 

10.45 

2260 

80 

218 

10.43 

2270 

223 

10.45 

2330 

90 

218 

10.42 

2270 

22? 

10.43 

2390 

100 

230 

10,40 

2390 

231 

10.42 

2410 

110 

238 

10.38 

2470 

229 

10,40 

2380 

120 

235 

10.36 

2430 

231 

10.38 

2400 
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TABLE  A.  3 (Continued) 
EXPERIMENTAL  STRAIN  AND  STRESS 

STRAIN  GAGE  NUMBER  S - 7 


TIME  CONFIGURATION  1 CONFIGURATION  2 


SEC. 

«•  106 

E-  iO-6 

(psi) 

e- 106 

E-IO'6 

(P*l) 

1G 

118 

10.5 

1230 

103 

10.5 

1130 

20 

155 

10.49 

163G 

150 

10.5 

1580 

30 

175 

» a 4-r 
1 v.  -rr 

1930 

i Oa 

10.48 

1900 

40 

205 

10.46 

2140 

199 

10.47 

2080 

50 

216 

10.45 

2260 

214 

10.46 

2240 

60 

218 

10.44 

2270 

2i5 

10.45 

2250 

70 

218 

10.43 

2270 

216 

10.43 

2260 

80 

219 

10.42 

2280 

224 

10.42 

2330 

90 

214 

10-  40 

2270 

223 

10.40 

/* 

4JIU 

100 

210 

10.37 

2180 

211 

10.38 

2190 

110 

209 

10. 35 

2160 

208 

10.36 

2150 

120 

205 

10. 32 

2120 

204 

10.34 

2110 
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EXPERIMENTAL  STRAIN  AND  STRESS 


STRAIN  GAGE  NUMBER  S - 8 

TiME  CONFIGURATION!  CONFIGURATION  2 


SEC. 

2<106 

E-10“* 

(ps!) 

e-lO6 

E-10"6 

(psl) 

10 

135 

10.5 

1420 

145 

10.5 

1520 

20 

209 

10.48 

2190 

2i! 

« A IA 

3U.  HT 

Win 

4.&4.V 

30 

258 

10,47 

2700 

258 

10. 48 

2700 

40 

291 

10.46 

3050 

291 

10.47 

3050 

50 

318 

10.45 

3320 

314 

10.46 

3280 

60 

334 

10.43 

3480 

328 

10.45 

3420 

70 

345 

10.41 

3590 

336 

10.43 

3510 

80 

354 

10.38 

3670 

348 

10.41 

3620 

90 

361 

10.36 

3740 

3 58 

10.40 

3720 

100 

380 

10. 34 

3930 

366 

10.37 

3810 

110 

390 

10. 32 

4020 

373 

10.35 

3860 

120 

398 

ia  on 

1 V.  ov 

4110 

380 

10.32 

3920 

!i  *1' ' 


TABLE  A.  3 (Contiruwi) 
EXPERIMENTAL  STRAIN  AND  SYRESS 

STRAIN  GAGE  NUMBER  S - 9 


TIME  CONFIGURATION  1 CONFIGURATION  2 


SEC. 

e.  106 

E-10"6 

(psl) 

e-106 

E-10"6 

<P«I) 

1C 

-135 

10,47 

-1410 

-150 

10.47 

-1570 

20 

-205 

10.45 

-2140 

“210 

10.45 

-2190 

30 

-*450 

10.42 

-2600 

-255 

10-43 

-2660 

40 

-279 

10.40 

-2900 

-2yj 

10. 40 

50 

-305 

10.36 

-3160 

-321 

10.37 

-3330 

60 

-318 

10. 33 

-3280 

-345 

10.34 

-3570 

70 

-326 

10,28 

-3350 

-358 

10.30 

-3680 

80 

-330 

10.24 

-3380 

-366 

10.27 

-3760 

90 

-331 

10.20 

-3380 

-368 

10.23 

-3760 

100 

-318 

10.16 

-3230 

-364 

10.18 

-3700 

T 10 

-307 

10.10 

-3100 

-353 

10. 15 

-3580 

120 

-288 

10.05 

-2900 

-339 

10.10 

-3420 

EXPERIMENTAL  strain  and  stress 


STRAIN  GAGl:  NUMBER  S - 11 

TIME  CONFIGURATION  1 CONFIGURATION  2 


SEC. 

e.  106 

E-10"6 

(p*l) 

e-106 

E-10"6 

(psl) 

10 

-269 

10.46 

-2810 

-254 

10.47 

-2660 

20 

-491 

10.42 

-5120 

— *t/  s; 

10.  43 

-4970 

30 

-639 

10.37 

-6620 

-639 

10.38 

-6630 

40 

-756 

10. 32 

-7800 

-768 

10.34 

-7940 

50 

-836 

10.26 

-8580 

-855 

10.28 

-8690 

60 

-886 

10.20 

-9040 

-921 

10.23 

-9420 

70 

-919 

10.15 

-9330 

-946 

10.17 

-9620 

80 

-933 

10.08 

-9400 

-969 

10.11 

-9790 

90 

-931 

10.02 

-9330 

-983 

10.04 

-9860 

<00 

"/*W 

A CO 
/•  /U 

-9270 

-983 

9.98 

-9810 

110 

-907 

9.87 

-8950 

-978 

9.90 

-9680 

120 

-678 

9.8Q 

-8610 

-961 

9.84 

-9460 
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